tDRmapper: challenges and solutions to mapping, naming, and quantifying tRNA-derived RNAs from human small RNA-sequencing data

Abstract Background Small RNA-sequencing has revealed the diversity and high abundance of small RNAs derived from tRNAs, referred to as tRNA-derived RNAs. However, at present, there is no standardized nomenclature and there are no methods for accurate annotation and quantification of these small RNAs. tRNA-derived RNAs have unique features that limit the utility of conventional alignment tools and quantification methods. Results We describe here the challenges of mapping, naming, and quantifying tRNA-derived RNAs and present a novel method that addresses them, called tDRmapper. We then use tDRmapper to perform a comparative analysis of tRNA-derived RNA profiles across different human cell types and diseases. We found that (1) tRNA-derived RNA profiles can differ dramatically across different cell types and disease states, (2) that positions and types of chemical modifications of tRNA-derived RNAs vary by cell type and disease, and (3) that entirely different tRNA-derived RNA species can be produced from the same parental tRNA depending on the cell type. Conclusion tDRmappernot only provides a standardized nomenclature and quantification scheme, but also includes graphical visualization that facilitates the discovery of novel tRNA and tRNA-derived RNA biology

Differential hepatitis C virus RNA target site selection and host factor activities of naturally occurring miR-122 3′ variants

In addition to suppressing cellular gene expression, certain miRNAs potently facilitate replication of specific positive-strand RNA viruses. miR-122, a pro-viral hepatitis C virus (HCV) host factor, binds and recruits Ago2 to tandem sites (S1 and S2) near the 5΄ end of the HCV genome, stabilizing it and promoting its synthesis. HCV target site selection follows canonical miRNA rules, but how non-templated 3΄ miR-122 modifications impact this unconventional miRNA action is unknown. High-throughput sequencing revealed that a 22 nt miRNA with 3΄G (‘22–3΄G’) comprised <63% of total miR-122 in human liver, whereas other variants (23–3΄A, 23–3΄U, 21–3΄U) represented 11–17%. All loaded equivalently into Ago2, and when tested individually functioned comparably in suppressing gene expression. In contrast, 23–3΄A and 23–3΄U were more active than 22–3΄G in stabilizing HCV RNA and promoting its replication, whereas 21–3΄U was almost completely inactive. This lack of 21–3΄U HCV host factor activity correlated with reduced recruitment of Ago2 to the HCV S1 site. Additional experiments demonstrated strong preference for guanosine at nt 22 of miR-122. Our findings reveal the importance of non-templated 3΄ miR-122 modifications to its HCV host factor activity, and identify unexpected differences in miRNA requirements for host gene suppression versus RNA virus replication

Model of fibrolamellar hepatocellular carcinomas reveals striking enrichment in cancer stem cells

The aetiology of human fibrolamellar hepatocellular carcinomas (hFL-HCCs), cancers occurring increasingly in children to young adults, is poorly understood. We present a transplantable tumour line, maintained in immune-compromised mice, and validate it as a bona fide model of hFL-HCCs by multiple methods. RNA-seq analysis confirms the presence of a fusion transcript (DNAJB1-PRKACA) characteristic of hFL-HCC tumours. The hFL-HCC tumour line is highly enriched for cancer stem cells as indicated by limited dilution tumourigenicity assays, spheroid formation and flow cytometry. Immunohistochemistry on the hFL-HCC model, with parallel studies on 27 primary hFL-HCC tumours, provides robust evidence for expression of endodermal stem cell traits. Transcriptomic analyses of the tumour line and of multiple, normal hepatic lineage stages reveal a gene signature for hFL-HCCs closely resembling that of biliary tree stem cells-newly discovered precursors for liver and pancreas. This model offers unprecedented opportunities to investigate mechanisms underlying hFL-HCCs pathogenesis and potential therapies

A Gene Expression Panel is Accurate for Diagnosis and Monitoring Treatment of Eosinophilic Esophagitis in Adults

Eosinophilic esophagitis (EoE) can be difficult to diagnose. We aimed to evaluate whether a gene expression score could differentiate adult EoE cases from non-EoE controls and to determine whether scores normalized after treatment for EoE

microRNA-122 Abundance in Hepatocellular Carcinoma and Non-Tumor Liver Tissue from Japanese Patients with Persistent HCV versus HBV Infection

Mechanisms of hepatic carcinogenesis in chronic hepatitis B and hepatitis C are incompletely defined but often assumed to be similar and related to immune-mediated inflammation. Despite this, several studies hint at differences in expression of miR-122, a liver-specific microRNA with tumor suppressor properties, in hepatocellular carcinoma (HCC) associated with hepatitis B virus (HBV) versus hepatitis C virus (HCV) infection. Differences in the expression of miR-122 in these cancers would be of interest, as miR-122 is an essential host factor for HCV but not HBV replication. To determine whether the abundance of miR-122 in cancer tissue is influenced by the nature of the underlying virus infection, we measured miR-122 by qRT-PCR in paired tumor and non-tumor tissues from cohorts of HBV- and HCV-infected Japanese patients. miR-122 abundance was significantly reduced from normal in HBV-associated HCC, but not in liver cancer associated with HCV infection. This difference was independent of the degree of differentiation of the liver cancer. Surprisingly, we also found significant differences in miR-122 expression in non-tumor tissue, with miR-122 abundance reduced from normal in HCV- but not HBV-infected liver. Similar differences were observed in HCV- vs. HBV-infected chimpanzees. Among HCV-infected Japanese subjects, reductions in miR-122 abundance in non-tumor tissue were associated with a single nucleotide polymorphism near the IL28B gene that predicts poor response to interferon-based therapy (TG vs. TT genotype at rs8099917), and correlated negatively with the abundance of multiple interferon-stimulated gene transcripts. Reduced levels of miR-122 in chronic hepatitis C thus appear to be associated with endogenous interferon responses to the virus, while differences in miR-122 expression in HCV- versus HBV-associated HCC likely reflect virus-specific mechanisms contributing to carcinogenesis. The continued expression of miR-122 in HCV-associated HCC may signify an important role for HCV replication late in the progression to malignancy

Transcriptomic Analysis of Chronic Hepatitis B and C and Liver Cancer Reveals MicroRNA-Mediated Control of Cholesterol Synthesis Programs

ABSTRACT Chronic hepatitis B (CHB), chronic hepatitis C (CHC), and associated hepatocellular carcinoma (HCC) are characterized by cholesterol imbalance and dyslipidemia; however, the key regulatory drivers of these phenotypes are incompletely understood. Using gene expression microarrays and high-throughput sequencing of small RNAs, we performed integrative analysis of microRNA (miRNA) and gene expression in nonmalignant and matched cancer tissue samples from human subjects with CHB or CHC and HCC. We also carried out follow-up functional studies of specific miRNAs in a cell-based system. These studies led to four major findings. First, pathways affecting cholesterol homeostasis were among the most significantly overrepresented among genes dysregulated in chronic viral hepatitis and especially in tumor tissue. Second, for each disease state, specific miRNA signatures that included miRNAs not previously associated with chronic viral hepatitis, such as miR-1307 in CHC, were identified. Notably, a few miRNAs, including miR-27 and miR-224, were components of the miRNA signatures of all four disease states: CHB, CHC, CHB-associated HCC, and CHC-associated HCC. Third, using a statistical simulation method (miRHub) applied to the gene expression data, we identified candidate master miRNA regulators of pathways controlling cholesterol homeostasis in chronic viral hepatitis and HCC, including miR-21, miR-27, and miR-33. Last, we validated in human hepatoma cells that both miR-21 and miR-27 significantly repress cholesterol synthesis and that miR-27 does so in part through regulation of the gene that codes for the rate-limiting enzyme 3-hydroxy-3-methyl-glutaryl-coenzyme A (HMG-CoA) reductase ( HMGCR ). IMPORTANCE Hepatitis B virus (HBV) and hepatitis C virus (HCV) are phylogenetically unrelated hepatotropic viruses that persistently infect hundreds of millions of people world-wide, often leading to chronic liver disease and hepatocellular carcinoma (HCC). Chronic hepatitis B (CHB), chronic hepatitis C (CHC), and associated HCC often lead to cholesterol imbalance and dyslipidemia. However, the regulatory mechanisms underlying the dysregulation of lipid pathways in these disease states are incompletely understood. MicroRNAs (miRNAs) have emerged as critical modulators of lipid homeostasis. Here we use a blend of genomic, molecular, and biochemical strategies to identify key miRNAs that drive the lipid phenotypes of chronic viral hepatitis and HCC. These findings provide a panoramic view of the miRNA landscape in chronic viral hepatitis, which could contribute to the development of novel and more-effective miRNA-based therapeutic strategies

Small tRNA-derived RNAs are increased and more abundant than microRNAs in chronic hepatitis B and C

Persistent infections with hepatitis B virus (HBV) or hepatitis C virus (HCV) account for the majority of cases of hepatic cirrhosis and hepatocellular carcinoma (HCC) worldwide. Small, non-coding RNAs play important roles in virus-host interactions. We used high throughput sequencing to conduct an unbiased profiling of small (14-40 nts) RNAs in liver from Japanese subjects with advanced hepatitis B or C and hepatocellular carcinoma (HCC). Small RNAs derived from tRNAs, specifically 30–35 nucleotide-long 5′ tRNA-halves (5′ tRHs), were abundant in non-malignant liver and significantly increased in humans and chimpanzees with chronic viral hepatitis. 5′ tRH abundance exceeded microRNA abundance in most infected non-cancerous tissues. In contrast, in matched cancer tissue, 5′ tRH abundance was reduced, and relative abundance of individual 5′ tRHs was altered. In hepatitis B-associated HCC, 5′ tRH abundance correlated with expression of the tRNA-cleaving ribonuclease, angiogenin. These results demonstrate that tRHs are the most abundant small RNAs in chronically infected liver and that their abundance is altered in liver cancer

Aminobisphosphonates reactivate the latent reservoir in people living with HIV-1

Antiretroviral therapy (ART) is not curative due to the existence of cellular reservoirs of latent HIV-1 that persist during therapy. Current research efforts to cure HIV-1 infection include “shock and kill” strategies to disrupt latency using small molecules or latency-reversing agents (LRAs) to induce expression of HIV-1 enabling cytotoxic immune cells to eliminate infected cells. The modest success of current LRAs urges the field to identify novel drugs with increased clinical efficacy. Aminobisphosphonates (N-BPs) that include pamidronate, zoledronate, or alendronate, are the first-line treatment of bone-related diseases including osteoporosis and bone malignancies. Here, we show the use of N-BPs as a novel class of LRA: we found in ex vivo assays using primary cells from ART-suppressed people living with HIV-1 that N-BPs induce HIV-1 from latency to levels that are comparable to the T cell activator phytohemagglutinin (PHA). RNA sequencing and mechanistic data suggested that reactivation may occur through activation of the activator protein 1 signaling pathway. Stored samples from a prior clinical trial aimed at analyzing the effect of alendronate on bone mineral density, provided further evidence of alendronate-mediated latency reversal and activation of immune effector cells. Decay of the reservoir measured by IPDA was however not detected. Our results demonstrate the novel use of N-BPs to reverse HIV-1 latency while inducing immune effector functions. This preliminary evidence merits further investigation in a controlled clinical setting possibly in combination with therapeutic vaccination

Model of fibrolamellar hepatocellular carcinomas reveals striking enrichment in cancer stem cells

The aetiology of human fibrolamellar hepatocellular carcinomas (hFL-HCCs), cancers occurring increasingly in children to young adults, is poorly understood. We present a transplantable tumour line, maintained in immune-compromised mice, and validate it as a bona fide model of hFL-HCCs by multiple methods. RNA-seq analysis confirms the presence of a fusion transcript (DNAJB1-PRKACA) characteristic of hFL-HCC tumours. The hFL-HCC tumour line is highly enriched for cancer stem cells as indicated by limited dilution tumourigenicity assays, spheroid formation and flow cytometry. Immunohistochemistry on the hFL-HCC model, with parallel studies on 27 primary hFL-HCC tumours, provides robust evidence for expression of endodermal stem cell traits. Transcriptomic analyses of the tumour line and of multiple, normal hepatic lineage stages reveal a gene signature for hFL-HCCs closely resembling that of biliary tree stem cells—newly discovered precursors for liver and pancreas. This model offers unprecedented opportunities to investigate mechanisms underlying hFL-HCCs pathogenesis and potential therapies

Combating subclonal evolution of resistant cancer phenotypes

Metastatic breast cancer remains challenging to treat, and most patients ultimately progress on therapy. This acquired drug resistance is largely due to drug-refractory sub-populations (subclones) within heterogeneous tumors. Here, we track the genetic and phenotypic subclonal evolution of four breast cancers through years of treatment to better understand how breast cancers become drug-resistant. Recurrently appearing post-chemotherapy mutations are rare. However, bulk and single-cell RNA sequencing reveal acquisition of malignant phenotypes after treatment, including enhanced mesenchymal and growth factor signaling, which may promote drug resistance, and decreased antigen presentation and TNF-α signaling, which may enable immune system avoidance. Some of these phenotypes pre-exist in pre-treatment subclones that become dominant after chemotherapy, indicating selection for resistance phenotypes. Post-chemotherapy cancer cells are effectively treated with drugs targeting acquired phenotypes. These findings highlight cancer's ability to evolve phenotypically and suggest a phenotype-targeted treatment strategy that adapts to cancer as it evolves